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ABSTRACT 


Hough  has  proposed  an  interesting  and  computationally  efficient  pro¬ 
cedure  for  detecting  lines  in  pictures.  In  this  paper  we  point  out  that 
the  use  of  angle-radius  rather  than  slope-intercept  parameters  simplifies 
the  computation  further.  We  also  show  how  the  method  can  be  used  for  more 
general  curve  fitting,  and  give  alternative  interpretations  that  explain 
the  source  of  its  efficiency. 
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Figure  Captions 


Fig.  1 
Fig.  2 
Fig.  3 


The  normal  parameters  for  a  line. 
Projection  of  colinear  points  onto  a  line. 
An  illustrative  example. 


I .  INTRODUCTION 


A  recurring  problem  in  computer  picture  processing  is  the  detec¬ 
tion  of  straight  lines  in  digitized  images.  In  the  simplest  case,  the 
picture  contains  a  number  of  discrete,  black  figure  points  lying  on  a 
white  background .  The  problem  is  to  detect  the  presence  of  groups  of 
colinear  or  almost  colinear  figure  points.  It  is  clear  that  the  problem 
can  be  solved  to  any  desired  degree  of  accuracy  by  testing  the  lines 
formed  by  all  pairs  of  points.  However,  the  computation  required  for  n 
points  is  approximately  proportional  to  n  ,  and  may  be  prohibitive  for 
large  n. 

Rosenfeld [1]  has  described  an^ingenlous  method  due  to  Hough [2]  for, 
replacing  the  original  problem  of  finding  colinear  points  by  a  mathemati¬ 
cally  equivalent  problem  of  finding  concurrent  lines .  This  method  involves 
transforming  each  of  the  figure  points  into  a  straight  line  in  a  parameter 
space .  The  parameter  space  is  defined  by  the  parametric  representation 
used  to  describe  lines  in  the  picture  plane.  Hough  chose  to  use  the 
familiar  slope-intercept  parameters ,  and  thus  his  parameter  space  was 
the  two-dimensional  slope-intercept  plane.  Unfortunately,  both  the  slope 
and  the  intercept  are  unbounded,  which  complicates  the  application  of  the 
technique .  In  this  note  we  suggest  an  alternative  parametrization  that 
eliminates  this  problem.  We  also  give  two  alternative  interpretations  of 
Hough's  method,  one  of  which  reveals  plainly  the  source  of  its  efficiency. 
Finally ,  we  show  how  the  method  can  be  extended  to  find  more  general 
classes  of  curves  in  pictures . 
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II.  FUNDAMENTALS 


The  set  of  all  straight  lines  in  the  picture  plane  constitutes  a 
two-parameter  family.  If  we  fix  a  parametrization  for  the  family,  then 
an  arbitrary  straight  line  can  be  represented  by  a  single  point  in  the 
parameter  space.  For  reasons  that  become  obvious,  we  prefer  the  so-called 
normal  parametrization.  As  illustrated  in  Fig.  1,  this  parametrization 
specifies  a  straight  line  by  the  angle  8  of  its  normal  and  its  algebraic 
distance  p  from  the  origin.  The  equation  of  a  line  corresponding  to  this 
geometry  is 

x  cos  0  +  y  sin  9  =  p 

If  we  restrict  0  to  the  interval  [0,tt),  then  the  normal  parameters  for  a 
line  are  unique.  With  this  restriction,  every  line  in  the  x-y  plane 
corresponds  to  a  unique  point  in  the  Q-p  plane. 

Suppose,  now,  that  we  have  some  set  {(x^y^,  ...,  (xn>yn)}  of  n 
figure  points  and  we  want  to  find  a  set  of  straight  lines  that  fit  them. 
We  transform  the  .points  (x^»jO  into  the  sinusoidal  curves  in  the  9-p 
plane  defined  by 

p  =  x.  cos  0  +  y.  sin  0  .  (1) 

i  i 

It  is  easy  to  show  that  the  curves  corresponding  to  colinear  figure 
points  have  a  common  point  of  intersection.  This  point  in  the  8-p  plane, 
say  (0o,po),  defines  the  line  passing  through  the  colinear  points.  Thus, 
the  problem  of  detecting  colinear  points  can  be  converted  to  the  problem 
of  finding  concurrent  curves . 
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FIGURE  1  THE  NORMAL  PARAMETERS  FOR  A  LINE 
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A  dual  property  of  the  point-to-curve  transformation  can  also  be 
established.  Suppose  we  have  a  setj^e^p  ),  <9k>Pk)}  of  points  in 

the  9-p  plane,  all  lying  on  the  curve 

P  =  Xq  cos  9  +  yQ  sin  9 

Then  it  is  easy  to  show  that  all  these  points  correspond  to  lines  in  the 
x-y  plane  passing  through  the  point  (x^jy^) .  We  can  summarize  these 
interesting  properties  of  the  point-to-curve  transformation  as  follows : 

1.  A  point  in  the  picture  plane  corresponds  to  a  sinusoidal 
curve  in  the  parameter  plane. 

2*  A  point  in  the  parameter  plane  corresponds  to  a  straight 
line  in  the  picture  plane. 

3.  Points  lying  on  the  same  straight  line  in  the  picture 
plane  correspond  to  curves  through  a  common  point  in  the 
parameter  plane . 

4.  Points  lying  on  the  same  curve  in  the  parameter  plane 
correspond  to  lines  through  the  same  point  in  the  picture 
plane . 

In  the  next  section  we  apply  these  results  to  the  problem 
of  detecting  colinear  points  in  the  picture  plane  and  show  how 
significant  computational  economies  can  be  realized  in  certain  situations. 

III.  APPLICATIONS  ANJ  ALTERNATIVE  INTERPRETATIONS 
Suppose  that  we  map  all  of  the  points  in  the  picture  plane  into  their 
corresponding  curves  in  the  parameter  plane.  In  general,  these  n  curves  will 
intersect  in  n(n  -  l)/2  points  corresponding  to  the  lines  between  all  pairs 
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of  figure  points .  Exactly  colinear  subsets  of  figure  points  can  be  found, 
at  least  in  principle,  by  finding  coincident  points  of  intersection  in 
the  parameter  plane.  Unfortunately,  this  approach  is  essentially  exhaus¬ 
tive,  and  the  computation  required  grows  quadratically  with  the  number  of 
picture  points . 

When  it  is  not  necessary  to  determine  the  lines  exactly,  the  computa¬ 
tional  burden  can  be  reduced  considerably.  Following  Hough's  basic  pro¬ 
posal,  we  specify  the  acceptable  error  in  9  and  p  and  quantize  the  Q-p 
plane  into  a  quadruled  grid.  This  quantization  can  be  confined  to  the 
region  0  ^  0  <  tt,  -  R  ^  p  <;  R,  where  R  is  the  size  of  the  retina,  since 
points  outside  this  rectangle  correspond  to  lines  in  the  picture  plane  that 
do  not  cross  the  retina.  The  quantized  region  is  treated  as  a  two-dimensional 
array  of  accumulators.  For  each  point  )  in  the  picture  plane,  the 

corresponding  curve  given  by  (1)  is  entered  in  the  array  by  incrementing 
the  count  in  each  cell  along  the  curve.  Thus,  a  given  cell  in  the  two- 
dimensional  accumulator  eventually  records  the  total  number  of  curves  passing 
through  it .  After  all  figure  points  have  been  treated,  the  array  is  inspected 

to  find  cells  having  high  counts.  If  the  count  in  a  given  cell  (9.,p.)  is  k, 

^  J 

then  precisely  k  figure  points  lie  (to  within  quantization  error)  along  the 

line  whose  normal  parameters  are  (9  ,p  ). 

i  j 

An  alternative  interpretation  of  the  point-curve  transformation 
can  be  obtained  by  recognizing  that  the  p  computed  by  (1), 

p  =  x.  cos  9  +  y.  sin  0  ,  (1) 

1  x 
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locates  the  projection  of  the  point  (x^jO  onto  a  line  through  the 
origin  with  slope  angle  9.  Thus,  if  a  number  of  figure  points  lie  close 
to  some  line  l ,  their  projections  onto  the  line  normal  to  1  are  nearly 
coincident  (see  Fig.  2) .  A  given  column  in  the  g-p  accumulator  array 
is  just  a  histogram  for  these  projections,  so  a  high  count  in  a 
given  cell  clearly  corresponds  to  a  nearly  colinear  subset  of  figure 
points.  A  variation  of  this  approach  was  used  by  Griffith [3]  to  find  long 
lines  in  a  picture. 

Let  us  investigate  how  the  computation  required  by  the  accumulator 
implementation  varies  with  the  number  of  figure  points.  To  be  more  speci¬ 
fic  about  the  quantization,  suppose  that  we  restrict  our  attention  to  d^ 
values  of  9  uniformly  spaced  in  the  interval  [0,tt)  .  Suppose  further  that 
the  p  axis  in  the  interval  [-R,  Rl  is  quantized  into  d  cells.  For  each 

a 

figure  point  (x^y^),  we  use  (1)  to  compute  the  d^  different  values  of  p 
corresponding  to  the  d^  possible  values  of  the  independent  variable  9. 
Since  there  are  n  figure  points,  we  need  to  carry  out  this  computation 
nd  times.  When  these  computations  are  complete,  the  d  d  cells  of  the 

-L  i  ^ 

two-dimensional  accumulator  are  inspected  to  find  high  counts.  Thus,  the 
computation  required  grows  linearly  with  the  number  of  figure  points . 
Clearly,  when  n  is  large  compared  to  d  ,  this  approach  is  preferable  to 
an  exhaustive  procedure  that  requires  considering  the  lines  between  all 
n(n  -  l)/2  pairs  of  figure  points. 
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FIGURE  2  PROJECTION  OF  COLINEAR  POINTS  ONTO  A  LINE 
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Another  alternative  interpretation  exposes  the  source  of  this 


efficiency.  Consider  again  Property  4  of  the  last  section:  Points 
lying  on  the  same  curve  in  the  9-p  plane  correspond  to  lines  through  the 
same  point  in  the  picture  plane.  When  the  curve  corresponding  to  figure 
point  (xi,yi>  is  "added"  to  the  accumulator,  we  are  really  computing  and 
recording  the  parameters  of  the  d^  lines  in  the  picture  plane  passing 
through  (x^,y^)  and,  because  9  is  quantized,  these  are  "all  the  lines  in 
the  plane"  passing  through  (x^,y_^)  .  Should  a  given  parameter  pair  ever 
recur  as  a  result  of  computing  the  d^  lines  through  some  other  figure  point, 
the  recurrence  will  be  reflected  in  an  increased  count  in  the  appropriate 
accumulator  cell.  Roughly  speaking,  then,  for  each  figure  point  the  quan¬ 
tized  transform  method  considers  only  the  set  of  all  d^  lines  through  that 
point,  whereas  more  exhaustive  methods  consider  all  (n  -  1)  lines  between 
the  given  point  and  all  other  figure  points . 

IV.  EXAMPLE 

The  following  example  illustrates  some  of  the  features  of  the 
transform  approach.  Fig.  3(a)  shows  a  television  monitor  view  of  a  box, 
and  Fig.  3(b)  shows  a  digitized  version  of  that  view.  A  simple  differ¬ 
encing  operation  locates  significant  intensity  changes  and  produces  the 
binary  picture  shown  in  Fig.  3(c).  This  120-by-120  picture  contains  many 
nearly  colinear  figure  points  that  can  be  fit  well  by  a  few  straight  lines. 

Sampling  0  at  d^  =  9  20°— increments  in  0  and,  quantizing  p  into 
d^  =  86  two-element  cells,  we  obtain  the  two-dimensional  accumulator  array 
shown  in  Table  I.  If  the  array  entry  at  (0Q, pQ)  is  kQ,  then  k^  figure 
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points  lie  on  parallel  lines  for  which  8  =0Q  and  p  lies  between  pQ  and 
Pq+2.  When  many  points  are  nearly  colinear*  the  entry  for  the  line  that 
fits  them  best  is  large.  The  largest  entry  in  the  table  occurs  at  (0°,-5) 

and  corresponds  to  the  middle  vertical  edge  of  the  box.  The  nine  circled 

entries  in  Table  I  correspond  to  locally  maximum  values  that  exceed  the 
arbitrary  threshold  of  35.  The  corresponding  nine  groups  of  nearly  colinear 
figure  points  are  shown  in  Fig.  3(d).  In  this  example,  it  happens  that  every 
group  corresponds  to  some  physically  meaningful  line  in  the  picture.  However, 
two  significant  lines  on  the  top  of  the  box  were  not  found,  one  because  it 
contained  very  few  points  and  the  other  because  it  fell  between  the  lines 

at  0  =  80°  and  9  =  100°.  The  20°  angular  quantization  interval  was  chosen 

to  keep  the  accumulator  array  small.  Clearly,  we  were  fortunate  to  have 
found  as  many  lines  as  we  did,  and  a  smaller  quantization  interval  -would 
:_have  to  be  used  in  practiiriT  ” 

A  few  remarks  concerning  some  limitations  of  the  transform  approach 
are  in  order.  First,  the  results  are  sensitive  to  the  quantization  of 
both  9  and  p.  Finer  quantization  gives  better  resolution,  but  increases 
the  computation  time  and  exposes  the  problem  of  clustering  entries  corres¬ 
ponding  to  nearly  colinear  points.  Second,  the  technique  finds  colinear 
points  without  regard  to  contiguity.  Thus,  the  position  of  a  best-fit  line 
can  be  distorted  by  the  presence  of  unrelated  figure  points  in  another  part 
of  the  picture.  A  related  problem  is  that  of  meaningless  groups  of  colinear 
points  being  detected.  In  our  example,  a  false  line  would  be  detected  if 
the  threshold  were  reduced  from  35  to  24,  the  value  needed  to  detect  the  top 
left-hand  edge  of  the  box. 
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TABLE  I  ACCUMULATOR  ARRAY  FOR  FIG.  3(c) 


Index  Terms 


Picture  processing,  pattern  recognition,  line  detection. 


An  important  special  use  of  the  transform  method  is  to  detect 
the  occurrence  of  figure  points  lying  on  a  straight  line  and  possessing 
some  specified  property.  For  example,  suppose  we  want  to  find  whether 
a  significant  number  of  figure  points  lie  on  a  line  through  the  point 
<x0,y0)  ^"n  Picture  plane.  As  we  have  seen  from  Property  4,  the 
normal  coordinates  of  any  such  line  must  lie  on  (or,  in  practice,  at 
least  near)  the  curve  p  =  xQ  cos  0  +  yQ  sin  0.  Hence,  the  transform 
process  can  be  carried  out  in  the  usual  way,  but  attention  can  be  re¬ 
stricted  to  the  region  of  the  0-p  plane  near  this  curve.  If  we  find  a 
cell  with  count  k  near  this  curve,  then  we  are  assured  that  k  figure 
points  lie  on  a  line  passing  (nearly)  through  the  point  (x0,yQ).  Sim¬ 
ilarly,  suppose  we  are  interested  only  in  lines  having  a  given  direction, 
say  0Q.  Again,  we  carry  out  the  process  in  the  usual  way,  but  restrict 
our  attention  to  a  subset  of  the  0-p  plane  in  the  vicinity  of  0  =  0Q. 

It  is  clear  that  the  general  transform  approach  can  be  extended  to 
curves  other  than  straight  lines.  For  example,  suppose  we  want  a  method 
to  detect  circular  configurations  of  figure  points.  We  can  choose  a  para¬ 
metric  representation  for  the  family  of  all  circles  (within  a  retina)  and 
transform  each  figure  point  in  the  obvious  way.  If,  as  a  parametric  rep¬ 
resentation,  we  describe  a  circle  in  the  picture  plane  by 

(x  -  a)2  +  (y  -  b)2  =  c2  , 

then  an  arbitrary  figure  point  (x^yj  will  be  transformed  into  a 
surface  in  the  a-b-c  parameter  space  defined  by 

(xi  -  a)  +  (y  -  b)  =  c 
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In  this  example,  then,  each  figure  point  will  be  transformed  into  a  right  cir¬ 
cular  cone  in  a  three-dimensional  parameter  space.  If  the  cones  corres¬ 
ponding  to  many  figure  points  intersect  at  a  single  point,  say  the  point 
(a0’b0’C0)’  then  fig^e  points  lie  on  the  circle  defined  by  those 

three  parameters.  As  in  the  preceding  case  of  straight  lines,  n©  saving  is 
effected  if  the  entire  process  is  performed  analytically.  However, 
the  process  can  be  implemented  efficiently  by  using  a  three-dimensional 
array  of  accumulators  representing  the  three-dimensional  parameter  space. 

In  principle,  then,  the  transform  method  extends  to  arbitrary  curves . 

We  need  only  pick  a  convenient  parametrization  for  the  family  of  curves  o^ 
interest  and  then  proceed  in  the  obvious  way.  A  parametrization  having 
bounded  parameters  is  obviously  preferable,  although  this  is  not  essential. 

It  is  much  more  important  to  have  a  small  number  of  parameters,  since  the 
accumulator  implementation  requires  quantization  of  the  entire  parameter 
space,  and  the  computation  grows  exponentially  with  the  number  of  parameters. 
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